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^ . We review the systematic uncertainties that have plagued attempts to ob- 

P^ . tain high precision and high accuracy from ground-based photometric measure- 

ments using CCDs. We identify two main challenges in breaking through the 
1% precision barrier: 1) fully characterizing atmospheric transmission, along the 
instrument's line of sight, and 2) properly identifying, measuring and remov- 



> 

in 



OO . ing instrumental artifacts. We discuss approximations and limitations inherent 

^ . in the present methodology, and we estimate their contributions to systematic 

O , photometric uncertainties. We propose an alternative conceptual scheme for the 

Q ! relative calibration of astronomical apparatus: the availability of calibrated de- 

r-| ! tectors whose relative spectral sensitivity is known to better than one part in 

'p! 10^ opens up the possibility of in situ relative throughput measurements, nor- 

P ! malized to a precision calibrated detector, using a stable but uncalibrated nar- 

c/3 ! rowband light source. An implementation scheme is outlined, which exploits the 

availability of tunable lasers to map out the relative wavelength response of an 
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imaging system, using a flatfield screen and a calibrated reference photodiode. 
The merits and limitations of this scheme are discussed. In tandem with care- 
5^ i ful measurements of atmospheric transmission, this approach could potentially 

lead to reliable ground-based photometry with fractional uncertainties below the 
percent level. 

Subject headings: atmospheric effects, instrumentation: photometers, instrumen- 
tation: detectors 



^Department of Physics and Harvard- Smithsonian Center for Astrophysics, Harvard University 
^Institute for Astronomy, University of Hawaii 



-2- 



1. Introduction and motivation 

It has proven difficult to achieve ground-based flux measurements with fractional uncer- 
tainties in accuracy of less than a few percent, using CCDs. This is much poorer than naive 
Poisson statistics would permit. Most photometry codes add a fractional flux uncertainty 
of a few percent in quadrature with the uncertainty from photon statistics, in recognition of 
various sources of systematic error (Stetson87), (Schechter93). The Sloan Digital Sky Survey 
(SDSS), for example, quotes (Stoughton et al 02) a systematic photometric zeropoint rms 
uncertainty of 2% across the survey flelds. Saha et al. (SahaOS) compare precision pho- 
tometry from multiple sources, and flnd systematic discrepancies at the 2-5% level. They 
attribute these discrepancies to passband differences and 

On the other hand, differential photometric measurements of sources with similar spec- 
tral energy distributions achieve (Hartman04) (Tonry05) precise differential flux measure- 
ments at the millimagnitude level, as long as the multiple sources can be captured in a single 
frame of the imaging system. Laboratory tests of the stability of CCD detectors indicate 
that they are capable of achieving Poisson-limited differential measurements at the level of 
one part in 10^ (Robinson95). The intrinsic stability of CCDs has led to ambitious plans for 
high resolution differential photometry from space (Borucki03). 

While very impressive relative flux measurements can now be made between objects that 
reside in a single exposure, determining the flux ratios between objects in different regions 
of the sky is much more difficult. This is one indicator that the temporal and directional 
variation in the optical transmission of the atmosphere are signiflcant limitations to precision 
photometry. 

It should be stressed that the normal procedures for photometric measurement are quite 
satisfactory for the purpose for which they were designed: measuring stars at the few percent 
level. Since 25% of stars are variable at the 1% rms level or more (and 50% at the 0.5% 
level) (Tonry05), and since stellar colors are often corrupted at the percent level by dust, 
there has been little impetus to do better. 

A number of forefront science issues demand better performance from ground-based sys- 
tems. The surge of interest in extrasolar planetary systems motivates searches for planetary 
transits, often via a subtle photometric signal. On the cosmological scale, using supernovae 
to probe the expansion history of the Universe requires making K-corrections (Peacock98) 
that rely upon detailed knowledge of the relative instrumental response as a function of 
wavelength. These considerations motivate taking a close look at how astronomical flux 
measurements are made, how the apparatus is calibrated, and how the data are analyzed. 
This paper will explore how we might break through the photometric precision barrier at 
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the percent level, for ground-based CCD flux measurements. 

In the sections that follow we will argue the following points: 

1. "Flatfielding" CCD images by dividing the flux detected in each pixel by a passband 
dependent sensitivity array is susceptible to systematic uncertainties due to differences 
in the spectral energy distributions of the flatfleld, the sky, and the source. 

2. Numerous approximations commonly made in the reduction of broadband images are 
likely to fail in attempts to measure flux with fractional uncertainty below the 1% level. 

3. It is possible and desirable to measure the relative throughput of the entire imaging 
system in situ, in which the wavelength-dependent throughput of the mirrors, corrector 
optics, fllter, and detector are determined using a tunable laser as a high brightness 
monochromatic source and a calibrated photodiode as a detector. 

4. There is considerable merit in obtaining ongoing measurements of atmospheric trans- 
mission over wavelength, using a dedicated spectrophotometric instrument in conjunc- 
tion with modern atmospheric models. 

5. Future imaging systems should be designed for ease of calibration and suppression of 
sources of systematic error. One obvious example is the inclusion of an independent 
determination of shutter timing. 

Numerous authors have explored the issues that pertain to the calibration of photomet- 
ric measurements, including (Bessell99), (Bohhn96), (Castelh94), (Colina94), (Fabregat96), 
(Fukugita96), (Hamuy92), (Hamuy94), (Hayes75), (Hayes75b), (Landoh92), (Magnier04), 
(Megessier95), (Oke83), (SahaOS), and others. Our focus here is less on the use of celestial 
calibrators than on exploiting the availability of calibrated silicon photodiodes, for charac- 
terizing astronomical apparatus, and on the diverse contributions to systematic uncertainty 
in astronomical flux measurements. 

Current practice in astronomical flux measurement makes an estimate of throughput 
for each passband, using fllter transmission data obtained from benchtop measurements, 
throughput estimates for optical components, and a detector quantum efficiency curve. 
Broadband measurements of standards are then used to "tweak" the system throughput 
function to force agreement between measurements and synthetic photometry (Bessell05). 
Unfortunately the integral measurement of flux over a broad passband does not uniquely 
determine what passband modiflcation is to be made. 
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The metrology chain that underpins most astronomical flux measurements dates back to 
an era when laboratory blackbody sources of molten metal were considered as fundamental 
standards. Heroic observing programs were undertaken to tie the spectrum of Vega to these 
laboratory sources. This accounts for the use of Vega as a "primary" spectral standard. In 
reality the primary metrology standard was the blackbody source on Earth, of course. 

The approach we suggest here exploits the availability of calibrated photodiodes as 
NIST-traceable metrology standards, with spectral response known at 10"^ level over the 
wavelengths relevant for CCD instruments. This allows us to calibrate the photon sensitivity 
of both imaging and dispersive apparatus as a function of wavelength, using a tunable 
narrowband light source. This technique allows us to measure relative system throughput 
for multiple passbands, thereby tying together the instrumental zeropoints across all filters. 
With this unambiguous system throughput data in hand, what remains is 1) to determine 
a single overall absolute throughput (or "effective aperture" ) for the system, common to all 
passbands and 2) ascertain the optical transmission properties of the atmosphere. 



2. The arithmetic of broadband photometry 

Most of the imaging measurements in astronomy use array detectors, with collecting 
optics and band-limiting optical filters, to measure the flux from sources of interest. Each 
pixel in the detector is illuminated by both astronomical sources and background radiation, 
termed "sky" in astronomical parlance. These different sources have different spectral energy 
distributions (SEDs). The signal Si from a source found in a pixel i is then a sum (over 
sources, plus sky) of integrals over wavelength. 



S^= Yl f ^jW X R.W X A dX, 



sources j 



where the J-'j are spectral photon distributions (SPDs, evaluated above the atmosphere) of the 
sky and all sources present, Ri{X) is the dimensionless system transmission of that specific 
pixel, including atmosphere, optics, filter, and detector, and Ai is the effective collecting 
aperture of the telescope for that pixel. (We adopt the convention here that the units of 
J-" are photons/nm/sec/cm^, derived from the conventional spectral energy distribution as 
J^ = F\ \/hc.) We call the product Ri{\) x A^ the "spectral aperture" of a pixel, it is that 
pixel's cross section for collection of light from a celestial object impinging at the top of the 
atmosphere. 

In general the transmissions Ri{\) are not constant from pixel to pixel, varying even as a 
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function of A (think of non-uniform interference filters, water spots on detector anti-refiection 
coatings, or the variations in internal interference giving rise to fringing). Conventional 
photometric analysis seeks to simplify this by renormalizing the broadband Si to a signal 
S with a known, common spectral aperture -R(A) x A, so that the resulting signal does not 
depend on which pixel i collected it. This is normally effected by division by a "flatfield", 
derived from a source of light which is intended to have no variation from pixel to pixel. It 
is obvious from Equation 1, however, that there is no single flatfielding correction, even in 
principle, that can be multiplicatively applied to the measured integral signal Si which can 
properly correct for the wavelength-dependent -Rj(A), for all possible SEDs. 

For example, an undulating supernova SED, a smooth galaxy SED, and the night sky 
emission lines might be all inextricably superposed on a single pixel. Because of internal 
fringing, a nearby pixel might have larger QE at a prominent sky emission feature, hence 
report a greater signal if the telescope were pointed to bring the feature onto this pixel 
(i.e. a fringe maximum instead of minimum). However, if the night sky aurora diminishes, 
this ratio between the two pixels would diminish, demonstrating that there is no unique 
flatfield scaling factor between these two pixels. Of course, sensible astronomers recognize 
the additive contamination of sky emission and subtract it, taking advantage of its temporal 
stability to expose a pixel to a common night sky SED but an ensemble average of celestial 
objects. The supernova and galaxy SEDs, are similarly disentangled by observing the galaxy 
when the supernova does not exist. While this illustrates important techniques of using time 
to improve differential photometry, it is important to bear in mind that in principle each 
pixel has its own sensitivity as a function of A. 

We argue below that it is possible to derive a common sensitivity function for an array 
of pixels which will be accurate at the sub-percent level for objects with SEDs which vary 
sufficiently slowly with wavelength (rapidly varying SEDs such as night sky emission can be 
handled by special procedures), but these considerations insist that fiuxes are best reported 
in such a "natural system" of the apparatus, as simply the observed integrated fiux over 
this common passband. Any SED-ignorant transformation to another photometric system 
through some set of linear equations (known as "color terms") is susceptible to systematic 
errors at the percent level or worse. With SED in hand, of course, it is simple to integrate it 
over the "natural system" and another passband to derive the signal which would have been 
collected in this passband. A clear example arises in supernova photometry, where using 
stars to derive transformation coefficients between photometric systems produces systematic 
errors due to differences in the SEDs of stars and supernovae. 

Saha et al. (SahaOS) and Sung & Bessell (SungOO) provide clear discussions of the 
importance of fully understanding the interplay between passband variations and SED vari- 
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ations. 



3. The arithmetic of hght collection and detection 

The light from a point celestial object arrives at the telescope nearly as a plane wave — 
the distribution of arriving photons is spatially uniform and essentially a delta function in 
angle. The telescope is a linear optical system which maps angles to positions on the focal 
plane, and so there is a transfer function H which describes how the intensity $ at positions 
X in the focal plane is related to incoming photon flux / at entrance aperture location x', 
pointing direction (3 of the optical axis, and angle o: relative to the optical axis: 



$(x, /3, t. A, P) = / /(x', a, /3, t. A, P) iJ(x, x', a. A, P) rfV d^oc. 



(2) 



The variables A and P represent wavelength and polarization (polarization is suppressed 
from now on). This integral expression can in principle be evaluated at any surface down 
the optical train. We flnd it most convenient to consider the plane just above the first 
mechanical or optical element, so that / from a point source in the sky is an unobscured 
plane wave of uniform surface intensity. 

It is approximately true (and a goal of optical design) that H be separable 

i7(x,x',Q;, A) Ri exp(iA;x-x7-F) exp(-iA;Q; • x') A(x') r7(A), (3) 

where k = 27r/A, F is the focal length (not to be confused with the notation for SED), A 
is the pupil transmission function (i.e. unity over the pupil, zero elsewhere), and 1]{X) is 
the wavelength response of optics and filter, but not including atmosphere or detector. The 
response $ at the focal plane is therefore the Fourier transform of the (potentially complex, 
phase-error carrying) input / at the entrance aperture, multiplied by a phase factor according 
to the direction of incidence. Ignoring the effects of diffraction, H will integrate to a non-zero 
response only when the argument of the exponentials is zero, 

X = Fa, (4) 

i.e. it maps angles to position, and picks up an amplitude factor of ~ cos a. 

Most contemporary astronomical detectors convert this incident photon fiux into pho- 
toelectrons where 

5(x)= /"<I>(x,A)g(x,A)rfA, (5) 

where Q is the detector's quantum efficiency, generally taken to depend only weakly on 
sub-pixel location, angle of arrival, or polarization. The response from a given pixel is the 
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area integral of S'(x) over the extent of that pixel, possibly convolved with diffusion within 
the detector or cross-talk from other pixels which are being addressed at the same time, 
run through an amplifier with a certain gain and linearity, and digitized by an A/D with 
its own behavior which may not be ideal at the 1% level. Understanding and suppressing 
imperfections in the signal chain is a fairly well defined electrical engineering problem, albeit 
one which merits careful consideration as we push for higher performance from the apparatus. 

The intensity which arrives at the entrance aperture is related to the spectral photon 
distribution T above the atmosphere 

/(A, a, /3, t) = T{\ ex, 13) T(A, a, (3, t), (6) 

where T is the atmospheric transmission, which depends upon wavelength, direction, and 
time. The challenge of ground-based photometry is to convert from the measured quantity, 
S'(x), namely the distribution of photoelectrons among the detector pixels, to ^(A), the 
photon spectra of the sources of interest. This requires properly understanding and correcting 
for 1) the transmission of the atmosphere T, 2) the optical transmission properties of the 
apparatus, H, and 3) the detector quantum efficiency Q. 

We advocate bundling the apparatus terms into the product HQ, which we should 
strive to measure in a way that minimizes systematic uncertainties, while considering the 
determination of the atmospheric transmission T as a distinct problem that merits its own 
dedicated monitor. We expect that the HQ product will not depend on where the telescope 
points (/3) and it should vary slowly with time, so it does not need to measured for every 
exposure. We also have full experimental access to the apparatus so we can measure HQ to 
arbitrary accuracy, at least in principle. The T term, however, should be measured for each 
exposure, with recognition that it bears some dramatic spectral features. 



4. Calibration of telescope and detector 

The availability of well characterized silicon photodiodes, with sensitivity vs. wavelength 
curves that can be traced to fundamental metrology standards (Larason98), provides us with 
the opportunity of exploiting these devices to perform a relative calibration of astronomical 
apparatus. This reference detector can be used in conjunction with a light source of arbitrary 
brightness to obtain a throughput vs. wavelength curve for astronomical apparatus. 

In principle, for each detector pixel we could create an appropriate monochromatic 
plane wave of known intensity at the entrance aperture of the telescope and thereby derive 
its spectral aperture function. This would also tell us the extent to which H does not 
integrate to a delta function in angle o: and detector position x, e.g. scattered or diffracted 



light, and it duplicates the illumination from celestial objects across the entrance aperture, 
so sensitivity in H to x' is removed. 

In practice, collection of such measurements would be too time consuming, so we endorse 
the normal procedure of filling the telescope entrance aperture with light from a fiat field 
screen and thereby simultaneously filling the cone of angles which illuminates the entire focal 
plane. However, we advocate using monochromatic light and monitoring the light level with 
a calibrated photodiode. Obtaining a series of flat fields at a discrete set of wavelengths and 
normalizing each fiatfield to the photon dose it received then provides us directly with a data 
cube which is an approximation of each pixel's spectral aperture. Another approach to this 
problem is described in (MarshallOS). 

This approximation fails to take account of breadth of (Fo: — x) ("scattered light"), and 
we do not have a better suggestion than to appeal to on-sky measurements to correct for 
this. We hope that the "scattered light" has the same spectral characteristics as the direct 
transmission part of H (i.e. if you paint your telescope structure red you will need to measure 
the scattered light in each of your filters), that it is spatially smooth, and temporally stable, 
so the fiatfield data cube can be multiplied by a smooth scaling function (which ideally is 
not a function of wavelength) to convert it to the accurate spectral aperture we seek. 

This "illumination correction" function is normally derived by rastering celestial objects 
across the detector during times that the atmosphere transmission is thought to be extremely 
uniform and constant. This correction is normally found to contribute at the few percent 
level, depending on the geometry of light scattering paths in the telescope and optics as well 
as the illumination pattern from the fiatfield screen. Manfroid (Manfroid95) (Manfroid96), 
and Magnier and Cuillandre (Magnier04) provide a thorough description of their experience 
in this regard. We obviously advocate a fiatfield illumination scheme which is extremely 
stable and uniform, and we next quantify how uniform it must be to achieve sub-percent 
photometry. 



4.1. Coupling between fiatfield illumination and piiotometry non-uniformities 

For the scheme outlined above, where a fiatfield screen is used to illuminate the entire 
focal plane, a number of effects can produce systematic errors in throughput measurements. 
We have made quantitative estimates of the corresponding requirements in illumination. 
Broadly speaking we can make errors in the spatial uniformity of the density of photons 
crossing the entrance aperture and we can make errors in their angular distribution. We 
take as given that an "illumination correction" will be performed; therefore many non- 
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uniformities which are not wavelength dependent will be calibrated out and not contribute 
systematic error. 

Observers ordinarily have to worry about matching the SED of the flatfield illumination 
to that of celestial objects (perhaps a 3000K incandescent filament with a color balance 
filter), but different portions of the flatfield screen may have significantly different SEDs (e.g. 
when illuminated by more than one flatfield lamp). By creating a data cube of normalized 
monochromatic flatfields we postpone this issue until we collapse the cube to a mean flatfield. 
This is discussed below. 

Spatial non-uniformities in illumination affect flatfield fidelity, but the ratio between 
photometry uniformity (i.e. flatfield fidelity) and spatial illumination uniformity is of or- 
der 5-10%, i.e. a relatively crude illumination density is tolerable. The extent to which 
"illumination correction" can mitigate these effects depends strongly on the resolution of 
the illumination correction, but we assume that this resolution is low. Effects which couple 
spatial non-uniformity to photometric non-uniformity include: 

• A pixel corresponding to an angle a receives light which is offset on the flatfield 
screen by cx.d, where d is the distance from the entrance aperture to the screen. The 
photometric uniformity is related to the integral illumination from these offset but 
partially overlapping patches on the screen. For the geometry of the MOSAIC imager 
at CTIO, for example, we estimate that we could tolerate large scale fractional surface 
brightness variation as large as SI /I ~ 15%. 

• The transfer function H may have spatially dependent terms (for example obscuration, 
dust on optics, degraded mirror coatings, or vignetting). Spatial non- uniformities in 
flatfield illumination thereby translate to spatial non-uniformities in detector illumi- 
nation. We estimate that the blurring effect of the integral over H reduce the latter 
to roughly 10% of the former, i.e. 10% illumination variations can still yield 1% pho- 
tometry, although this depends on details of the nature of the obscuration and where 
it is. 

• Our scheme depends on a flux normalization from a photodiode which does not exactly 
monitor the light seen by a typical pixel. We are then susceptible to normalization 
errors if this mismatch is wavelength dependent. 

The angular distribution within the entrance cone which maps to the detector translates 
directly to photometric accuracy, although the "illumination correction" procedure alluded 
to above can in principle correct any non-uniform angular pattern which is stable. This 
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distribution must therefore be uniform (or known) at the sub-percent level for sub-percent 
photometry. 

• Spatial non-uniformities of the flatfield screen angular response are another case of 
spatial non- uniformity described above. Flatfield screens are often lit at a glancing 
angle, and the scattering function between the illumination direction and the (nearly) 
normal direction into the telescope can vary markedly at extreme angles. 



• 



Spatial non-uniformities across the screen can also cause bandpass errors. We normally 
use interference filters in a converging beam which undergo a wavelength shift which 
depends on the angle of incidence: 

A(^) = Ao(l - sin^ e/nlf/^ ~ Ao(l - ^'/S)- (7) 

For example, illumination of the CTIO 4-m flatfield screen by a Lambertian source 
at the top of the prime focus has a radial intensity fall-off across the screen which 
depopulates the outermost cone angles. At 700 nm, the average wavelength shift would 
be 1.2 nm, as compared with 1.5 nm for a uniformly filled cone, so such a flatfield passes 
through a filter which is 0.3 nm redder than that encountered by a celestial object. 

• A "thousand points of light" flatfield screen consisting of 1000 downward-pointing fibers 
distributed randomly across the entrance aperture is not only spatially non-uniform, 
but would require very careful injection of light since fibers "remember" the angular 
distribution of light fed into them for a great distance. 

The angular distribution outside the detector entrance cone will create scattered light 
which may differ from that from the sky. The desire for a very uniform angular distribution 
usually leads to an extremely wide illumination pattern, generally the full 27r steradians. 
Effects which affect the flatfield include: 

• The angular distribution of scattered light from celestial objects is usually limited by 
the dome aperture, whereas a flatfleld screen radiates at all angles, giving rise to quite 
different scattering contributions. Attempts to mitigage this effect by limiting the 
spatial size of the screen to the entrance aperture have very little utility in a modern 
dome which closely flts the telescope. The "illumination correction" is essential for 
flatfleld screens which radiate signiflcantly outside of the entrance cone. 

• Observers are all familiar with the dramatic smears of scattered light occuring when 
moonlight illuminates the telescope structure. Similar effects occur from sources of light 
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in the dome other than the flatfield screen, for example dayhght leaking in through 
cracks or instrumentation LEDs. We argue below that it is vital that flatfield illumi- 
nation be shutterable, so that a "fiatfield" is the difference between the lit and dark 
screen. This differential technique will be effective to the extent that temporal variation 
in light leaks in the dome are slow compared to the on-off sequencing time. Cloudless 
days are therefore optimal for obtaining these calibration data. 

In summary we desire the angular uniformity of our flatfield source to match the pho- 
tometric accuracy desired (or at least stable at that level), but the spatial uniformity can be 
~ lOx less uniform because obscuration couples in only at the 5-10% level or so. Scattered 
light is likely to be a signiflcant contributor no matter what we do, implying that application 
of an "illumination correction" function is required. We therefore want to try to maxmize 
the smoothness of the scattered light contribution and minimize its wavelength dependence. 



4.2. Implementing this Calibration Scheme 

In order to characterize the system throughput we need to provide three components: 
a source of monochromatic light, a means of injecting it into the telescope aperture, and a 
way to measure its intensity. 



4-2.1. Monochromatic light source 

An 0.25" pixel requires some 10^ photons for a flatfleld, but subtends only 2 x 10~^^ of 
27r steradians, so at 20% efficiency a light source must provide 0.1 joule to a 2tt steradian 
injection system. While a monochromator can do this, the tradeoff between spectral purity 
(we need ~1 nm to see fringing) and exposure time (acquisition of a 500 wavelength data 
cube through N filters is some lOOOA^ exposures) is very challenging unless we can deliver 
the light into much less than 2it steradians. 

Another scheme is a pulsed tunable laser (the Vibrant by Opotek is one example) which 
offers many advantages. These devices can deliver tens to hundreds of mW of power into an 
optical fiber, tunable across the wavelengths of interest. The short coherence length of the 
light avoids any speckle effects. We have installed such a unit at CTIO and are pleased with 
the results (Stubbs05). 
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4-2.2. The Flatfield Screen 

As outlined above, we wish to provide a full-pupil illumination source that has surface 
brightness uniform to 10%. Most flat field systems are reflective: one or more sources 
project light towards a screen that then reflects it back into the telescope system. A simple 
reimaging system with a single on-axis source will exhibit a radial intensity dependence 
of roughly cos^{6) where 6 is the angle from the telescope axis to the screen element dA, 
measured at the source. Minimizing this requires a large standoff distance, in order to reduce 
the angle 9. 

Modern telescope enclosures are engineered so as to minimize the enclosed volume, and 
the top end of the telescope is often precariously close to the inner wall of the enclosure. 
This lack of stand-off distance makes it very difficult to achieve unobscured, uniform surface 
brightness illumination with a reflective scheme. 

We are in the process of developing a uniform back-lit flatfield screen. An optical flber 
is placed between a reflective backing and a diffusing screen. The optical flber is engineered 
to produce substantial light leakage along its length, and the layout of the flber in the 
sandwich is designed to produce uniform surface brightness over the surface of the screen, 
compensating for the variation in brightness along the flber's length. Tests we have carried 
out on prototypes of this approach are very promising, and we intend to implement this 
scheme for both the Pan-Starrs and LSST systems. 



4-2.3. The illumination monitor 

In order to map out the spectral response of the telescope we must have a means of 
measuring the light intensity from the flatfleld screen. This can be readily achieved with a 
calibrated photodiode. 

Ideally the photodiode would measure the light emanating from the entire extent of 
the flatfleld screen within the same solid angle destined for the focal plane. This may 
not practical, however, so the photodiode's view of the flatfleld may not sample the entire 
spatial extent, and/or it may see a larger solid angle than does the focal plane. Obviously 
an effective design depends on the spatial and angular uniformity of the flatfleld emission, 
so we advocate analysis and experimentation to flnd the best compromise for photodiode 
illumination, given the dominant systematic errors of a given flatfleld implementation. We 
have had good success in using the photodiode in a simple pinhole camera conflguration, 
monitoring the light emanating from the entire surface of the flatfleld screen. 
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4- 2.4- A focal plane monitor? 

We note that placing another cahbrated photodiode in the focal plane is valuable for 
separating the throughput of the telescope from the QE of the detector. Aging aluminum 
or filters can thereby be distinguished from CCDs onto which outgassed substances are 
condensing. To the extent that the CCD QE is stable (which is likely to be the case over 
timescales of months), this also provides a measure of the gain of the CCD signal chain. 



5. Measuring the optical properties of the atmosphere 

A trio of phenomena dominate the transmission properties of the atmosphere (SlaterSO) 
(Houghton77). The physical processes ^ that determine the atmospheric transmission T(A) 
include: 

1. Rayleigh scattering from the atoms and molecules of atmospheric gases, 

2. Small particle scattering from aerosols suspended in the atmosphere, and 

3. Molecular absorption, in particular by O2, O3, and H2O. 

In addition to the imperfect transmissive properties of the Earth's atmosphere, at wave- 
lengths longer than 7500 A the atmosphere is a significant source of narrowband emission, 
mostly from OH molecular transitions in the ionosphere. 

The atmospheric transmission can be a complicated function of airmass, particularly 
at wavelengths where molecular absorption is a major factor. These absorption features 
are responsible for the "telluric" absorption lines in ground-based spectra, and as they ap- 
proach saturation the airmass dependence is a complicated function of airmass (Adelman96), 
(Houghton77), (Guiterrez82), (SchusterOl). Assuming broadband transmission with a simple 
scaling with airmass is inappropriate for precision compensation for atmospheric absorption, 
a fact that is fully appreciated among infrared astronomers (Manduca79). 

With the atmospheric transmission factored out as a separate measurement problem, we 
can use a dedicated piece of apparatus to monitor celestial sources in order to measure T(A). 
This can be done with multiband photometry, with a dispersive instrument, or both. We 



^Scintillation can also produce photometric variation for systems with aperture D comparable to the 
atmospheric Fried parameter tq. We will assume all integrations are long enough to suppress this effect to 
below the 1% level. 



-14- 



propose making measurements of atmospheric transmission at the time, and along the same 
hne of sight, as the primary instrument is being used. This will require separate stand-alone 
apparatus. 

The Pan-Starrs project intends to use a combination of a wide-field single-chip imager 
to monitor transparency across the entire FOV of the imager, in conjunction with a spec- 
trograph to monitor the wavelength dependence of atmospheric throughput. The wide-field 
imaging aspect has been in routine use at CFHT for some time (Cuillandre02). Here we will 
focus on the dispersive measurement. 

Fortunately, much is known about the optical properties of the atmosphere (SlaterSO), 
so extracting a parametric description of its transmission qualities should not be difficult. 
Unfortunately we know surprisingly little about the angular and temporal correlations of 
atmospheric emission and transmission, so the conservative approach is to ensure that the 
T system is keeping pace with the main telescope. 

One advantage to using a dispersive instrument to monitor atmospheric transmission 
is that the emission spectrum of the sky is also then available. Since many astronomical 
sources of interest are significantly fainter than the sky, in should in principle (knowing the 
system response -R(A)) be possible to estimate the number of sky counts each pixel i receives, 
and explicitly subtract this from the sum in equation (1). 

We intend to use a modest telescope and a commercial fiber-fed low dispersion spectro- 
graph to monitor appropriate sources within or close to the fields being observed with the 
Blanco 4m. We will explore how to best extract a parametric description of the 3 processes 
of concern, namely Rayleigh scattering, aerosol scattering, and molecular absorption. We 
intend to use these parameters as input to a comprehensive computer model of the atmo- 
sphere (such as Modtran), in order to allow us to correct for molecular absorption features 
that are well below our spectral resolution. 

Initial progress along these lines is described by (Granett05). 



6. Image Processing Implications and Opportunities 

6.1. Mean spectral response and fiatfield construction 

Ignoring issues such as sensitivity to polarization and photon arrival direction, the con- 
struction of a monochromatic fiatfield from our fiatfield datacube is simple in principle, al- 
though the removal of the scattered light component which might differ in intensity or SED 
from celestial light is an experimental challenge. We normally do not observe in monochro- 
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matic light, however, but through relatively broad filters with AA/A ~ 0.2, and we have 
argued above that it is therefore impossible to correct the fiux ratios reported by different 
pixels to a consistent answer for all celestial sources. On the other hand, the mathematically 
consistent approach of reporting the observed pixel responses along with a data cube of 
response functions for each pixel (effectively the illumination-corrected, normalized flatfield 
data cube) is not satisfactory either. We do recognize that there are celestial sources such as 
planetary nebulae or HII regions whose emission is essentially a delta function in wavelength 
for which this is the only approach which will work. (Even then caution is necessary: we have 
alluded to the wavelength shift of a fiatfield arising from non-uniform screen illumination 
that could cause large error for an emission line which falls right at the edge of a bandpass.) 

One need only glance at the image of a monochromatic planewave to appreciate that 
there is indeed a very substantial variability in the pixel-to-pixel response. Fringing in the 
near IR causes a response oscillation in Q every d\ = \^ /{2nw) or equivalently a change in 
thickness by \/{2n), where w is the thickness of the detector and n ~ 3.5 is the index of 
refraction of silicon. Diffraction around dust particles also creates a wavelength-dependent 
response of H at all wavelengths. These effects are reduced by a fast f/ratio and a broad 
cone of incoming light, but they are usually present at the few percent level at least, and a 
thin device with an AR coating not optimized for the IR can fringe at the 10% level or more. 

The response to celestial sources is the integral over wavelength of the pixel response 
function Ri times the celestial SPD. If we wish to use this to infer what a system with a 
mean response R would be to this SPD, we will make a fractional error which is 

6S_ ^ J J^jX) {R{X) - R.{X)) dX 
S J J^{X) Ri{X) dX • ^ ^ 

Obviously the error will be large for objects where changes in ^(A) coincide with deviations 
Ri{X)—R{X), integrated over the bandpass. Conversely, if we can keep i?j(A)— -R(A) relatively 
small, ^(A) relatively smooth, and the bandpass is relatively wide, the error will be small. 
These three factors can be traded off against one another: if we were willing to keep track 
of each pixel's Ri{X) there would be no error; trying to go from a pixel's Ri to some grossly 
different bandpass Rgtd will incur larger error; objects with smooth J^{X) will have more 
accurate results than those which vary, etc. Integration of equation (8) for "reasonable" 
SPD and detector response over a "reasonably" wide bandpass indeed yields errors which 
are less than 1 percent. 

We consider a "reasonably" smooth SPD to have no more than isolated 20% variation 
per nm (e.g. H&K break or molecular band edges), and otherwise 5% variation per nm 
(e.g. supernovae). A "reasonable" ratio Ri/ R between pixel response and mean response 
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does not exceed a 2%/nm derivative at bandpass edge, a slope of 2% across a bandpass, and 
oscillations of 10% with period 2.5 nm in regions redder than 750 nm, where sky emission 
lines become important. A "reasonably" wide bandpass is AA/A ~ 0.2. We regard objects 
which do not meet this criterion to be "photometric challenges" and analysis seeking sub- 
percent photometry should pay attention to the individual Ri lest an unfortunate coincidence 
of emission or band edge with response oscillation cause large photometric error. We regard 
pixels which do not meet this criterion to be "photometrically challenged" and should be 
flagged as such. We assert that it is sufficient, for most apphcations of interest to astronomers, 
to construct a single system response function -R(A) per passband, along with a single 2-d 
flatfield image that corrects for pixel-to-pixel variations and illumination nonuniformities, 
appropriately weighted across the passband. Both of these can be extracted from the stack 
of monochromatic flatfield images. 

We advocate applying the multi-pronged technique used by most astronomers: subtrac- 
tion removal of sources with violently varying SPDs (emission lines), division by the fiatfield 
data cube which has been integrated over a mean bandpass or a fiatfield obtained from a 
continuum source, reporting a fiux value whose collection pixel signature has been suppressed 
to below the 1% level, and supplying a mean response function vs. wavelength for the entire 
array. 

The image processing stages that we consider appropriate, using the approach described 
here, include: 

• Begin by correcting the raw fiatfield data cube, wavelength by wavelength: 

— Remove ambient light contamination by subtracting source on/off pairs. 

— Normalize each fiat to a uniform photon dose by dividing by the integrated sig- 
nal seen by the calibrated photodiode, after accounting for its photon detection 
efficiency. Correction for a spatially non-uniform shutter response is appropriate 
here. 

— Divide out large-scale unwanted artifacts, such as edge dimming from the Jacobian 
of the distortion (we want to preserve fiuxes after division by fiatfield), but leaving 
in desired features such as edge dimming from vignetting. 

— Calculate the illumination correction, by integrating a nominal bandpass fiatfield 
(described below), dividing it into images of a star field taken under photometric 
conditions at many different offsets, measuring the fiuxes and using all ratios for 
each star to estimate a model (e.g. polynomial or coarse grid) for the illumination 
correction. 
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— Multiply all the normalized wavelength slices by the illumination correction model. 

• Construct a mean response function R{X) for all pixels by averaging the responses 
Ri{X) of individual pixels over the array. (Each Ri should be scaled to the same value 
first, we are not interested in wavelength independent dimming from effects such as 
vignetting or distortion.) 

• Verify that differences between the mean response and that of any pixel, Ri/ R — 1, 
does not exceed a 2%/nm derivative at bandpass edge, a slope of 2% across a bandpass, 
and oscillations of 10% with period 2.5 nm in regions redder than 750 nm, where sky 
emission lines become important. 

• Verify that the integration of source SPDs of interest, over a bandpass of width 100 nm, 
have location-dependent errors relative to the mean bandpass of < 1%. Detectors 
which are worse or SPDs which are more variable will have larger non-uniformity. (It 
is possible to construct counterexamples of allowed SPDs and responses which will not 
meet the resulting 1% criterion.) 

• There is a degree of freedom in how we construct our 2-d flatfield, from the normalized 
monochromatic flats. Celestial objects approximate blackbodies between 1000 K and 
infinite temperature, hence have SPDs (photon rate per nm) which go approximately 
as A+^° (Wien exponential) to A^'' (4500 K) to A^^ (Rayleigh- Jeans tail) across the 
optical region of the spectrum. Across a broad filter bandpass of AA/A ~ 0.2 this 
amounts to a difference of 40%. More relevant, our requirements on the constancy 
of Ri/ R — 1 will give us 1% photometry for any of these SPDs, but it is important 
to recognize that our photometric accuracy is in the instrumental system with mean 
response function i?(A). 

• To the extent the flatfield matches the celestial SPD, the pixel-to-pixel variation will be 
removed perfectly. We therefore should select some intermediate weighting function, 
e.g. none (A"'"'') or filament temperature 2800 K (A+^), and integrate the flatfield 
datacube with that weight. This is then reflected in our mean, dimensionless response 
function. Another option is to weight the monochromatic flats so as to correspond to 
a particular calibration spectrum of interest, such as AB magnitudes (A~^), although 
that doesn't automatically confer additional accuracy in AB magnitudes. 

• It is acceptable to illuminate the flatfield screen with white light (and thereby do 
an analog integration of the flatfield datacube) but the division of image by flatfield 
impresses the SPD of the flatfield on the data, so the flatfield SPD is a component of 
the system response R and must be measured. Apart from the advantages of measuring 



• 
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the system spectral response, we regard measurement of a white light SPD and keeping 
it stable to be difficult enough that we believe that it is preferable to accumulate and 
integrate the monochromatic data cube. 

If we were to estimate a flux in a "standard" fllter which differs from our actual 
bandpass without any correction for the source's SPD we could clearly make large 
errors to the extent that the bandpasses differ. However, the traditional method of 
using "color terms" to estimate an ideal response from an observed one is clearly an 
approximation based on a linear (or sometimes non-linear) flt to a one-parameter family 
of SPDs, usually stars or blackbodies. These cannot be used for 1% photometry of 
SPDs which differ by more than 10% from the SPD which generated them. E.g. a 
set of color terms created from stellar SPDs are adequate for 1% photometry of most 
stars, possibly for galaxies, but not for supernovae. 

It is possible to do a posteriori improvement of photometry if we can estimate the SPD 
of an object of interest, for example, a supernova at a given age and redshift. We can 
then return to the datacube of pixel response functions and integrate the pixel response 
and mean response against both the supernova SPD and the chosen flatfleld SPD, and 
use the ratios to correct the flux from the flattened image to that which would have 
been reported by pixels which follow the mean response exactly. 



6.2. Explicit Sky Subtraction? 

One could imagine determining the sky emission spectrum using the data from a ded- 
icated dispersive sky monitor and a model for OH emission. The contribution of the sky 
spectrum to the counts in each pixel could then be subtracted using the monochromatic 
flatfleld stack. The fringe pattern that arises from strong sky lines illuminating a detector 
with a spatial variation in sensitivity is normally very stable in character but can be quite 
variable in amplitude. It remains to be seen whether the sky emission data will provide a 
fringing correction that is superior to that obtained from stacking images and subtracting 
an appropriately scaled fringe frame, where the amplitude of fringing in each image must to 
be determined before the master fringe frame can be scaled and subtracted. 



6.3. Photometric measurements 

We can hope to have reasonably accurate absolute photometry if we know that the sky 
has a well behaved transparency, for example if the conditions are "photometric" as judged 
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by the sky monitor. However, photometry at the 1% level is unlikely to be possible without 
recourse to celestial photometric standards. The present situation of a handful of calibrated 
photometric standard stars is currently improving via SDSS to photometrically calibrated 
stars found in every field of view. 

The advent of the astrometric-photometric survey of Pan-STARRS which should begin 
in 2006 will improve the situation yet more. This will extend photometric measurements to 
the 1% level in a very well characterized set of filters to the entire sky visible from Hawaii. The 
accurate characterization of the Pan-STARRS bandpasses and knowledge of the standard 
object's SPDs should permit use of these object's calibration to bring any observation to an 
accurate value for physical flux integrated over the bandpass in use. 

As detectors become larger and photometry requirements become more demanding, it 
becomes necessary to regard the sky transparency as a variable across the field of view. This 
can potentially be gauged from the sky monitor, but a more likely procedure will be to use 
the very high density of calibrated objects to map the variation in transparency. 



7. Considerations for the design of new apparatus 

Many aspects which limit the precision of photometry are legacies of the days before 
electronics, computers, and large scale detectors, such as the use of a single "exposure time" 
instead of a spatially variable shutter timing function. Some improvements can be made 
with a few new calibration procedures and subroutines, such as linearity corrections. Others 
require new apparatus. 



7.1. Signal chain electronics: accuracy and stability 

No digital data will be better than the analog electronics which produces it. A CCD 
output MOSFET dissipates about 10 mW which is enough power to cause local warming. 
A large signal can change this power by 10%, and a CCD amplifier can be sensitive enough 
to temperature that this changes its gain. This would show up as non-linearity (which 
might be different for a bright star than a flatfield), but of course this non-linearity could be 
sensitive to the overall temperature of the CCD. The near-IR QE of CCDs is very sensitive 
to temperature as well, as much as 1%/K at 1 /im. This strongly suggests that CCD 
temperatures be regulated to no worse than 1 K. 

CCD amplifiers are typically non-linear near the 1% level, and the full well of the CCD 
pixels is often not reached because the amplifier has become so non-linear that the data are 
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useless. There are many satisfactory techniques to measure the hnearity of a CCD amphfier 
and we recommend that this (a) be done on a regular basis, (b) a correction be routinely 
applied to the data stream, and (c) the results be available as meta-data. If the linearity 
changes with time it is a sign that that amplifier and signal chain should not be trusted for 
accurate photometry and it is important to know this. 

The gain of an amplifier and signal chain is a related parameter. In principle it is 
removed by flatfield division, but in practice the gain of an ailing CCD amplifier or signal 
chain A/D can also change rapidly. There are also techniques to monitor this gain at the 
few percent level and again, it should be done regularly and available as meta-data so that 
a misbehaving channel is flagged before it is used to do high precision photometry. 



7.2. Filter Passband Width 

In order to ensure that pixel-to-pixel variations in Ri{\) do not produce spurious results, 
it is important that the optical passbands be wide enough to encompass many periods of 
sinusoidal response variation with wavelength. The source SED plays a role here, of course. 
In general we estimate that passbands of AA/A ~ 0.2 should be sufficient to keep position- 
dependent flux errors below 1% for most sources of interest. 



7.3. Shutter timing 

Many observers pay attention to the fact that shutter timing is not accurately reported 
by the instrument, and that shutter timing depends on position in the focal plane; others do 
not. Mechanical shutters typically have non-uniformities at the 20 msec level, which trans- 
lates to a 1% photometric error for a 2 sec exposure. Of course this is a common exposure for 
flatfields and bright photometric standards. It is straightforward to determine empirically 
the shutter timing function over the detector by comparing short and long exposures if the 
shutter is reproducible, but if the shutter's performance depends on orientation or temper- 
ature this will not work. We strongly advocate that (a) shutters be equipped with sensors 
to report the shadow trajectory across the detector to an accuracy of 0.1 msec, and (b) this 
trajectory be coded into an exposure's metadata and reduction pipelines be equipped to use 
a "shutter timing function" for reductions rather than a single exposure time. 

A recent example of shutter timing issues apparently producing systematic discrepancies 
in photometry is described in (StetsonOS). 
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7.4. Designing for calibration ease 

Pan-Starrs is an ambitious system nearing deployment; LSST and SNAP are in the 
design phase. For many applications of the data that these systems will produce, systematics 
are likely to limit science. This implies that the calibration aspects should be afforded a 
high priority, including: (1) building in a good uniform flat field screen, properly baffled, 
(2) designing for fast readout so that acquisition of flatfield data is efficient, (3) allowing for 
more sophisticated flatfielding and calibration processing in the reduction pipelines. 



8. Conclusions and Next Steps 

A different approach to the challenge of ground-based photometry, taking advantage of 
well-calibrated detectors as metrology standards, could well allow measurements with relative 
flux uncertainties well below the one percent level. An important ingredient is measuring 
directly the optical properties of the atmosphere. 

It is our intention to pursue these ideas as part of the ESSENCE supernova survey, 
as well as in the Pan-Starrs and LSST projects. We are currently building and deploying 
equipment for the CTIO Blanco 4-m telescope and 8k mosaic as part of the ESSENCE pro- 
gram. We also expect to implement most of these ideas on the flrst Pan-STARRS telescope, 
scheduled for flrst light in mid-2006. 
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